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Abstract 

We analyse Ricci flow (normalised/un-normalised) of product manifolds -unwarped as well as 
warped, through a study of generic examples. First, we investigate such flows for the unwarped 
scenario with manifolds of the type x S™, S" x H™ , H™ x H" and also, similar multiple products. 
We are able to single out generic features such as singularity formation, isotropisation at particular 
values of the flow parameter and evolution characteristics. Subsequently, motivated by warped 
braneworlds and extra dimensions, we look at Ricci flows of warped spacetimes. Here, we are able 
to flnd analytic solutions for a special case by variable separation. For others we numerically solve 
the equations and draw certain useful inferences about the evolution of the warp factor, the scalar 
curvature as well the occurrence of singularities at flnite values of the flow parameter. We also 
investigate the dependence of the singularities of the flow on the initial conditions. We expect our 
results to be useful in any physical/mathematical context where such product manifolds may arise. 
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I. INTRODUCTION AND OVERVIEW 

Ricci and other geometric flows jlj] have been a topic of active research interest in both 
mathematics and physics, over the last few years. Though introduced many years ago by 
Hamilton in mathematics and in the work of Friedan sl in the context of a-models and 
string theory, it was Perelman's proof of the Poincare conjecture using Ricci flows jj] which 
generated most of the current interest. Much of the subsequent work has thus focused on 
using ideas associated with Ricci flows, in the context of general relativity and gravitation 
jsl. In this article, we follow a similar line of thought with reference to metrics on warped 
manifolds of various types. 

The definition of a warped manifold in mathematics is somewhat different from that 
used in physics. In mathematics, product manifolds with different overall conformal factors 
dependent only on a certain parameter, are termed as warped. However, in physics, warping 
refers to an extra coordinate dependence of one or more of these conformal factors-in other 
words, the conformal factor associated with one manifold in the product when assumed to 
be dependent on one or more of the coordinates of the other manifold, is said to be warped. 
Otherwise, it is termed as unwarped. We shall follow this latter definition in our work here. 

We restrict our attention to Ricci flow of specific types of product manifolds. Earlier work 
on such manifolds can be found in . We begin, following the work on non-Bianchi types in 
the penultimate section in 0], by looking at products of the form S" x S'", §" x H™, x H" 
for arbitrary m and n, and with different conformal factors (dependent on a parameter A) 
associated with the canonical metric on each manifold in the product. This study essentially 
involves unwarped products in the sense described in the previous paragraph. We also look 
at multiple products (more than two) as a possible generalisation. What do we look for 
in our investigations? Apart from exact solutions of the flow equations (which may not be 
possible always), we focus attention on the following aspects: (a) blow- up/dropping to zero 
value of any of the conformal factors at a finite A, (b) isotropisation (i.e. equal value of the 
conformal factor at a finite A) (c) fixed points/curves of the flows (d) dependence of flow on 
initial conditions. 

In the second half of our paper, we look at specific warped product manifolds largely 
inspired by currently fashionable ideas in the physics of extra dimensions . In particular, our 
product manifolds here are essentially similar to those used extensively in warped braneworld 
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models introduced in the work of Randall and Sundrum, about a decade ago 8j. A Ricci 
flow of such manifolds seems to provide interesting pointers in the physics of braneworlds-an 
aspect we discuss towards the end of our article. 

We have studied normalised Ricci flows in the unwarped cases while for the warped 
ones we primarily investigate the unnormalised Ricci flow. In the former, we are mainly 
concerned with Riemannian metrics, whereas in the latter we discuss the Lorentzian case in 
detail because it is the one which is physically relevant. 

A word here about the essential mathematical motivation of Ricci flows may not be too 
inappropriate. In any method to find the 'best' metric on a manifold the first thing to do 
it to define a natural evolution of metrics. One such evolution of metrics w.r.t. a 'time' 
parameter (not the physical time, in any sense) is the Ricci flow. Once we have defined the 
flow it is important to prove that the flow exists for all time and converges to a geometric 
limit. In a case the flow does not converge, the corresponding metrics degenerate and one 
then needs to relate the degeneration with the topology of the manifold. Among many good 
candidates, a specific choice is an evolution equation of vector fields (rather, a family of 
vector fields) on the space of metrics known as the Ricci flow equation and defined as: 



dgij 



-2% (1) 



dX 

For the normalised Ricci flow, the corresponding equation turns out to be: 

^ = -2i?., + ^(i?)^7., (2) 

where (R) = "j^r- Thus the normalised flow ceases to be different from the unnormalised 
one if we consider non-compact, infinite volume manifolds (with a finite value for J RdV). 
Moreover, for constant curvature manifolds the second term in the normalised flow equation 
reduces to -R. In this article, we shall first look at normalised flows in the next section. 
Subsequently, we shall deal with un-normalised flows. 



II. RICCI FLOW OF UNWARPED PRODUCT MANIFOLDS 



As mentioned in the Introduction, in this Section, we focus attention on unwarped product 
manifolds. We discuss each case separately by studying the evolution equations and solving 
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each dynamical system analytically. We present our results through figures which illustrate 
the nature of evolution of each conformal factor in the corresponding product manifold. In 
this Section, our study is based on normalised Ricci fiows, since all our manifolds are of 
constant curvature. 



A. On S"(/D) X S'^(/E) 

Here we focus on product manifolds of the aforementioned type with D and E representing 
the conformal factors for the and in the product. D and E are functions of the fiow 
parameter A. 

Let us compute the various curvature quantities(Ricci tensor, Ricci scalar) of this product 
manifold. We denote the metric on S"'(y^) by ag^^ where g^^ is the canonical metric on 
n -sphere. In the same way, the metric on the product manifold S"(V^) x S™(V^) looks 
like: g = Dg^n + Egs^,. Let Y and V be the unit vector fields on S" and S™ respectively 
and X be a unit vector field either on §" or which is perpendicular to both Y and V. It 
can be shown that (by Koszul's formula) (VyX, V) = as [Y, X] is either zero or tangent 
to S" and likewise with [X, V]. Thus VyX = if X is tangent to S™. Also VyX is tangent 
to §" if X is tangent to S", showing that VyX can be computed on S"(-\/U). This shows 
that if X, r are tangent to §" and U, V are tangent to §'^(/E), then yK{X AV) = where 
m is the curvature operator. Also 9\{X AY) = ^X AY and 9\{U AV) = AV. In 
particular, all sectional curvatures lie in the interval [0,max{-^, -^}]. From this we can say: 
Rc{X) = 

Rc(V) = ^^~^^ V or we can write Rcg = {n — l)(7§n + (m — l)5'§m and R = "'-^^^ + 
Hiil^zll^ Therefore we can make inference that S"(-\/zJ) x §™(-\/^) always has constant 
scalar curvature and it is an Einstein manifold exactly when = -^^^ (this is possible 
either n, m > 2 ot n = m = 1) . It is also worthwhile to mention that this product manifold 
has constant sectional curvature only when n = m = 1. Now we have all the ingredients in 
our hand to write normalised Ricci fiow equations. These are given as: 

dD , , 2 {n[n — \) m(m — 1) 1 , , 

d\ ^ ^ (n + m) \ E ] ^ ' 

^ = -2(m-l)+ , ' | M^-l) ^ m(m-l) j^ 
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The above two equations represent a first order dynamical system. For different n and m 
we get different dynamical systems. For n,m = 1 we have D and E as constant and the 
metric is a flat metric on §^ x S^. 



1. General Solution 



After some simplification, we have: 

dD _ 2m{n - 1) 2m(m - 1) D 
dX m + n {n + m) E 



(5) 



dE _ 2n{m - 1) 2n{n - I) E 
d\ m + n [n + m) D 

We may analyse the following general form for the above pair of equations. Similar 
equations will arise for the other product manifolds to be discussed later on. Solutions can 
be found for arbitrary a, 6, c and d. Specific values can be inserted for each case separately. 
Thus, we now analyse the generic system: 

dx , X dy V 

_ = -a + 6 - ; — = -c + d - (7) 
dX y dX X 

Dividing one equation with the other, we get: 

dx y{yd — cx) 



dy x{hx — ay) 
which can be solved by putting y = vx: 



(8) 



dv v'^ia + d) — vic + b) v^ci + VC2 ,„s 

X — — ^ (9) 

dx {h — av) {h — av) 

Defining ci = {d + a) and C2 = — (c + 6), the general solution of x and v will look like: 

— ln( ) ln(cit' + C2) = Inx — In & (10) 

C2 CiV + C2 Cl 

For §"(v//^) X §™(/E), with Ci = 2{n - 1) and = -2(m - 1) we have: 

Tn 

InxH Inw = In fc, ^ x^y"" = (11) 

m + n 
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If we put the above back in the original equations, we get the profile of x or y. 

dx b 

This can be integrated to get A(x) in terms of hypergeometric functions. For n = m or 
n = l,m, as we shall see one can obtain x and y as functions of A. Otherwise, in general, 
the functions A(x) or X{y) are not invertible. We now discuss a few illustrative, special cases. 



2. Special cases 

a. S\VD) X S^i^/E) 

Here n = 1 and m = 2. The equations become: 



dD _AD dE _ 2 
^ ~ 3 E ' dA ~ ~3 



(13) 



The solutions are straightforward: 



^0 Y ^..^ ^ 2 



^^^^ = ^^\-E;^) ' E{\) = Eo--\ (14) 

where at A = 0, -D = -Dq, E = Eq. One notices that D increases with increasing A and 
diverges to infinity at a finite A = |i?o- At the same value of A, E drops to zero. The flow 
therefore ends at this value of A-a curvature singularity appearing because of the shrinking 



of the part to zero radius. Fig 1(a) shows the aforementioned behaviour. 



In general, for products of the type S^(v-D) x S''"(v£') we end up with the following 
expressions for -D(A) and E{\). 

(\ m 
^o-fe^A; (m+l) 
The conclusions for arbitrary m, are similar to those for m = 2. 
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(a)L»o = Eo = l for x §' 
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(b)Do = 2, £^0 = 3 for x S' 
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(c)i:>o = 1, Eo = 0.9999 for x §2 (d)i:)o ^ 0.9999, = 1 for x 

FIG. 1: D{\) (dashed) and E{X) (continuous). The horizontal axes is A. 



b. §2(^) X sWe) 



Here, we turn to the more interesting case of E>'^{\/D) x S^(V^). The equations are: 



dD _ D dE _ E 

'dX~~ ^E ' dA~~ 



(16) 



We easily note that DE = constant = c{say). Putting -| = -f the first equation looks like: 



dD _D^-c 
dX c 



(17) 



and the solution with, D{X = 0) = Dq is: 



D{X) = 



£'o(l + exp^) + v/^(l - exp^) 

2A 

LDo(1 -exp^) + v^(l +exp^) 



(18) 
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From DE = c one can find E{X) which is given as: 



^(A) = 



Dnil 



exp 



+ yc(l + expv^; 



Doil + exp^-) + v/c(l - exp 



2A 



(19) 



Thus, it is clear that Eq = For Dq = Eq = ^c, we get D{X) and -E(A) as constant for 
all A. For Dq > ^ we have a divergence in -D(A) while -E'(A) drops to zero at the same 
A = \d- On the other hand, if Dq < -\/c, the behaviour of -D(A) and E{X) are reversed. In 
either case the value of A at which the divergence occurs is given as: 



A. = ^lnf^^ 



(20) 



Fig, 1(c) and 1(d) demonstrate these features. 



It may be noted that the above results for S^(V^) x S^(V^) holds in the more general 
setting, for §"'(-\/U) x S"(V^) with an appropriate scaling of A. The general equations are: 



dD 



(n 



DM 4) 



dE 



[n 



(21) 



The solution for -D(A) turns out to be: 



D{X) 



2(n-l)A 2(7i-l)A ^ 

-Do(l + exp ^ ) + -yc(l — exp ^ ' 



2(n-l)A 2(n-l)A 

Do(l-exp ^ ) + 75(1+ exp ) 



(22) 



Further A^ is given as: 



Ad = — 7Tln(- 



(23) 



2(n-l)""^Do- 

It is worthwhile to note that in this case, D = E =constant represents a fixed point 
(soliton) of the flow. 



3. Numerical Evaluations 

Though we have exact solutions formally, it is easier to analyse the nature of the flow if we 
obtain numerical solutions of the dynamical systems involved. This is because the analytic 
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solutions is the cases other than the ones quoted above are necessarily non-invertible. 

We shall focus here on cases where n ^ m. Let us see if it is possible to have -D(Ao) = 
E{\q) for n 7^ m. Looking at the general equations one can evaluate the derivatives assuming 
this {D = E) to be true. It turns out that, 

'dD_ 

Thus for n = m there is no chance of D to be equal to E at some A. On the other hand, this 
is always possible for n ^ m. We illustrate this fact by solving the equations numerically, 
as a dynamical system. 



An 



2m(m — n) 
n + m 



dE 



. An 



2n{n — m) 
n -\~ m 



(24) 



a. X §3(/E) 



The equations in this case are: 



dD _ _6 12D dE _ _8 AE 

'dX ~ ~5E ' ~dX~~5^5D ^ ^ 



Fig 1(b) below demonstrates the facts mentioned above. 

Similar results may be obtained for any other case with n ^ m. In summary, we observe 
the following: 

• There are fixed curves defined by = -^^—4. 

• As one of the radii (D or E) increases (or becomes singular), the other drops to zero. This 
will also be apparent from alnD vslnE plot, which will follow the equation: 

InD = -^\nE + ^\nk. 

n n 

• Reversing the initial conditions can lead to a reversal of the behaviour of D or E. This is 
easily seen in the case n = m where the equations go over to each other if we swap D and 
E. However, this is also true, as we have checked (not shown here) for n ^ m. 

• Generically similar behaviour for all m, is also seen when one of the two manifolds is 
and the other is (with m 7^ 1). 
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B. On §"(7(1?)) X I]r{y^{E)) 

Let us now turn towards analysing another type of product manifold (as mentioned 
above). For this case, the metric looks like: g = Dg^n + Eg^-m where gs^ is the canonical 
metric on n dimensional sphere and (^e™ is the metric on a m dimensional hyperbolic space. 
The normalised Ricci flow equations turn out to be: 

^ = -2(n-l)+, ' {n{n-l) _m{m-l)^ 

dE , , 2 [n[n — \\ m(m — 1) 1 . 
— = 2 m - 1 + <^ ^ ^ ' ^— '- \ E (27) 

where, as before, D and E are functions of the flow parameter A. 



1. General Solution of this system 

As in the previous subsection, after some simplification we have: 

2m{n — 1) 2m(m — 1) x . 2n{m — 1) 2n(n — 1) y , . 

m + n [n + m) y m + n [n + m) x 

Following the general analysis in terms of x, y and the a, b, c, d given earlier, we note that 
the values of a, b, c, d, for this case are: 

(n — 1) , (m — 1) (m — 1) , (n — 1) , , 

a = 2mi { ; b = -2m) { ; c = -2nj { ; d = 2ni (29) 

(m + n) (m + n) [m + n) (m + n) 

Also, Ci = {d + a) = 2{n — 1), C2 = — (c + 6) = 2{m — 1). Notice the difference in signs as 
compared to the S" x case. Further, the relation between x and y will be: 

TTl 

lnx + \nv = Ink ^ x"" y'^ = k (30) 

m + n 

Though this appears to be the same as for x §"^, a difference will appear in the explicit 
solution of x{\) and y(A). For instance, 

dX b m + n / rt-t \ 
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and, consequently, one can find y{X). Let us look at some special cases now. 



2. S^{^/D) xM^{VE) 

In this rather simple case, the difference quoted above is easily seen. The equations are: 

dX 3E ' dX 3 ^ ^ 



The solution is: 



^^^^=^^\^;;fjx) ' ^(^) = ^o+3A (33) 



where at A = 0, D = Dq, E = Eq. Note the change in sign, which suggests that there is no 
singularity in the evolution of D, starting from A = 0, unlike the previous case of x S^. 
The result here can easily be generalised for §^ x 



rrm 



3. S"(V£>) X W^y/E) 

Here, the equations turn out to be: 

dD , , ^ X -D dE , / - N /„.\ 

_._(„_1)_(„_1)_ ; _ = („_1) + („_1)_ (34) 

The features to note here are: 

• There are no fixed curves, unlike the previous case for S'" x . 

• There is always a possibility of having D = E a.t some A. 

• As before DE is a constant, but evolution of -D(A) and E{X) differs as compared to the 
previous cases. 

C. M" X M"* 

Finally we turn to products of the form H" x H"*. Here, the metric takes the form 
g = Dg^u + Eg^m where gm^{m) is the canonical metric on n(m) dimensional hyperbolic 
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space. The normalised Ricci flow equations become: 

dD , , 2 fnfn — 1) m(m — 1) "1 , , 

— = 2 n - 1 <^ ^ ^ ' + '- \ D 35) 

dE , , 2 { nin — 1) m(m — 1) 1 „ , , 

— = 2 m - 1 <^ ^ ^ ' + \ E 36) 

dX ^ ' n + m\ D E } ^ ' 

The analysis for these equations follows in the same way as for the earlier cases. We list 
below are few characteristic observations. 

• The flxed curves are given as = -^^-this is similar to the x S'" case. 

• The possibility oi D = E ai some A exists as long as n 7^ m. 

• The flow ends at singular points, depending on initial conditions. 



D. Multiple products and general trends 

What could happen if we take multiple products of the or H" or mixed products of 
various kinds? In general, it is difficult question to address. However, we can surely try out 
some specific examples. We illustrate things using one such example, below. 



1. Sp{VF) X S"(/D) X S'"(/B) 

For this triple product, the metric looks like: g = Fg^p + Dg^n + Eg^m where g§n is 
the canonical metric on n dimensional sphere. The normalised Ricci fiow equations are as 
follows: 

dF , , 2 fpfp— 1) nin — V] m(m — 1) "1 , , 

^ = -2(^-1) + - ? [V{v-r) n{n-X) m{r.-l)\^ 

d\ ^ ' {j>^n^m)\ F D E J ^ ' 

dE , , 2 fnfn— 1) nin — V) m(m — 1) 1 

— = -2m-l + - r <^ i-^^-— ^ + ^ + -\E 39 

d\ ^ ^ + n + m) \ F D E } ^ ^ 

The above fiow equations can be analysed in a way similar to the cases discussed earlier. The 
dynamical system is far more complicated. A few general comments can be made without 
solving the equations. 
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(a)Fo = 3, - 4, So = 5 (b)Fo ^3,Do^2,Eo^3 

FIG. 2: F(A) (dashed), D{X) (thick) and E{X) (thin) for x x S'^. The horizontal axes is A 

• Ifn = m = p= l then F, D and E do not evolve. 

• If any one of n,m or p is equal to one (say p = 1), then the system of equations for D, E 
can be solved independently as was done for the simple products discussed in the previous 
subsections. These solutions can be used to find the evolution of F. 

• li n = m = p then the equations are simpler, but unlike the case of a product of two 
manifolds, we do not have any fixed curves here. lin^m^poTn = my^p there are fixed 
curves which can be found by just putting the R. H. S. of the above equations to zero and 
solving the algebraic equations. 

• An intersection of all three functions D, E and F at some A is possible as long as 
p < n < m, n"^ + m? > p{n + m) and p^ + < n{j) + m). This can be checked by using 
D = E = F on the R. H. S. of the general equations and assuming that ^ > 0, ^ < at 
the same A. One can easily show that the above restrictions are impossible to achieve. This 
is also visible in the numerical example given below. 



2. S2(/p) X §3(/D) X S^{Ve) 

In this case the equations are: 

dF _ 14 4F 8F dD _ 8 AD 8D 

~dX~~Y^3D^3E ' ~d\ ~ 9F^ 3E ^ ^ 

dE _ 10 4 E 4 E 

IX ~ ~Y ^ 9D^ 3F ^ ^ 

We numerically solve this system. Fig|2] shows some typical evolution features for F, D and 
E. 
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The figures illustrate the fact that there can be a pair among D, E and F which can 
intersect at some A. Different pairs intersect at different points. Also, there can also be a 
situation without any intersection, this seems to be the case when a pair among D. E and F 
have the same initial value. Furthermore, divergences appear and the flow ends at a finite 
A, as was observed in the earlier cases with a product of two manifolds. 

It is possible to generalise the above results to products of arbitrary number of manifolds 
of differing dimensions. One can easily write down the general evolution equations by inspec- 
tion. The case with products of S*" and and also only can also discussed. Generic 
features are possibly not very different. Furthermore, instead of taking the normalisation 
on the whole manifold, one may consider, following related work in [3], cases with volume 
normalisation along sections. We will discuss some of these issues in a later article. For 
the moment, we prefer to switch over and move on to the perhaps more interesting case of 
warped products. 

III. RICCI FLOW OF WARPED PRODUCT MANIFOLDS 

In spacetime models with extra dimensions, braneworld scenarios provide many ex- 
amples of product manifolds. We consider one such case where the manifold is a product 
of Minkowski spacetime M. and the real line TZ^. In addition, we have warping in the sense 
described earlier in the Introduction, which makes the ensuing analysis different from the 
ones described in the previous sections. 

In order to consider the evolution under Ricci flow of the warped product Ai x TZ^ we 
assume the line element of the form: 

ds^ = e2/(^'^) {~dt^ + dx^ + dy^ + dz^) + rl{a, X)da^ (42) 

Note that the functions / and are assumed to be functions of a as well as a non-coordinate 
parameter A. We now evolve the metric functions in the above line element according to a 
un-normalised Ricci flow, with A as the flow parameter. 
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Th relevant geometric quantities for this metric are - 

g2/ / /■yx 
- Rtt = R.. = Ryy = R,, = ^ /" + 4/'2 - (43) 

R^^ = -Aif" + f' 

R=-±(^2f" + 5r 
After some straightforward algebra, the Ricci flow equations become - 

f" + 4/'2 _ f<\ (46) 
rc J 

f" + f" - ^) (47) 

here " / " represents |^ and " /' " represents ^ and similarly for Tc. The equation for 
/ may be thought of as a certain type of nonlinear heat equation with the equation for r^ 
appearing as a constraint. 

We make note of the remarkable fact that, the above Ricci flow equations are independent 
the signature of the part of the metric that contains a scaling with e^-'^. 

We now try to understand the behaviour of / (especially as a function of A) by solving the 
flow equations. The resulting solutions will tell us how the metric, as well other geometric 
quantities, evolve along the flow. 

A. Separable metric functions: exact results 

Let us first assume that 

r,(cr,A) = r,(A) ; /(a. A) = /.(a) + /^(A) (48) 
This assumption enables us to write the EqJl6] and EqHT] in a variable separable form as: 

rlfx = (/; + ^f'J) = K, (49) 



(44) 
(45) 
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(50) 



For the cr-part we get 



which imply Ki = —K2 and - 



K2 
12 



(51) 



(52) 



Therefore, we can easily solve the A evolution to get, 



(53) 



where, we make use of Ki = —K2 and the initial condition that at A = we have / = and 
Tc = 1 that is no scaling in the metric. 

To simplify, we can further write Ki = l/2Ac and thus the final metric turns out to be: 



a; 



exp 



/2X 



-dt^ + dx"^ + dy'^ + dz^) + da"^ 



(54) 



The Ricci scalar for the abovestated metric evolves as, 



R 



5/2 



A + Ac 



(55) 



We can prove that the solution obtained above is the only solution of the set Eq|l6] and 
Eqj47l when Tc = rc(A). Under these assumptions Eql47l can be written as: 



4 (/" + /'2) = 4B{Xf 



(56) 



This readily yields = 8 f B{\ydX = jS; and then solving for /, 



/i = ±5(A)a + C(A) 
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or 

/2 = A{X) + In [cosh {B{X)a + C{X))] 
Using / = /i in EqHHl then gives - 

(3 X (^±Ba + = 4B^ 

This imphes B = constant and thus / = ±Ba + C(A). This is just the variable separable 
case considered earlier. 

Similarly, for / = /2 we have - 



P 



A + tanh {Ba + C) x (jBa + 6^ = [l + 3 tanh^ {Ba + C)] 



Since, both sides represent functions analytic in a, we compare their Taylor series around 
the point cr = 0. The comparison of first 5 Taylor coefficients gives - 

B = -—sech^C ; C = —- tanh C (l + tanh^ C) 

jj jJ 

-B^ + A(3 + 3B^ tanh^ C = 

AB^ tanh Csech'^C = 

(4 - 3 cosh 2C) sech^C = 

The only solution to these is i? = i.e. / = y4(A) + ln [coshC(A)] and again / is separable. 
This is a limiting case of the earlier solution when Ac oo and leads to a static flat metric. 
Thus, EqinHis the general (and only) solution of the Ricci flow equations, when = rc{\). 



B. Nonseparable metric functions: numerical results 

To deal with a more general situation, we now drop the assumption, Tc = Tc(A) and 
consider the full non-separable case of EqsJlHl&llZl We observe that if /(cr. A) and rc{a, A) 
are a solution of the flow then so are f{aa, a^A) and rc{aa, a^A). Thus, we look for scaling 
solutions of the form f{y) and rc{y) with y = Note that it is also possible to work 
with - instead of y, though, with such a choice, large y will correspond to small A and vice 



17 



versa. We further assume that either cr > or ctq < cr < (Ji, i.e. we consider, as the extra 
dimensional space the open half line or an open interval. 

In terms of the new variable y and defining u = ^, the flow equations reduce to the 
following first order dynamical system, 



Also, from Eq|l5]we can write - 



d?f 2 {df\ 5 (dfV 1 fdf\ fdr, 



dy"^ y \dy J 2 \dy J \dy J \ dy 



-— + 12y^A\^ 

y 



(60) 



We shall use this latter expression in our evaluations. 

Note that the above dynamical system is invariant under scaling transformations of the 
type {y ay, f ^ f,A ^ A/a,rc arc). We can, thus, choose the scale by specifying 
initial conditions at y = 1 and solutions for conditions at other values of y can be obtained 
simply by the appropriate scaling. Hence, starting from y = 1 we can study the evolution 
of the system in the future {y > 1). 

Further, as all our results will be in terms of y = ^ we must therefore fix a (at some 
o" 7^ 0) in order to understand the evolution in A at that a and subsequently repeat things 
for another a. Physical considerations dictate that Vc > and since the equations are 
independent of / we can take /(I) = 1. We first choose the initial conditions at y = 1 as 



{f,A,rc) = (1,±1,1). The results of the numerical integration are plotted in Fig, 3(a) 



3(d) 



We note that when ^4(1) = +1, / eventually goes to a constant and drops to zero in 
the future {y > 1). Thus the flow exists for arbitrarily large values of y i.e. in the future and 



18 



f(l)=l,A(l)==l,r(l)==l 



f(l)=l, A(l)==-l,r(l)==l 



f(l)=l,A(l)==l,r(l)==l 



1.000 : 
0.995 : 
0.990 : 
0.985 : 
0.980 : 
0.975 : 
0.970 : 
0.965 - 



1 



(a)/ v/s y 




1 i A 

y 



f(l)=l, A(l)==-l,r(l)==l 



f(l)=l,A(l)==l,r(l)==l 



{h)A v/s y 



f(l)=l,A(l)==-l,r(l)==l 




{c)rc v/s 2/ 



f(l)=l, A(l)==l,r(l)==l 



f(l)=l, A(l)==-l,r(l)==l 




{i)(7^R v/s y 



FIG. 3: Results of numerical integration of dynamical system. Figures on left are for A[\) = +1 
and those on right are for A{\) = —1 
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the geometry tends towards a singular geometry only as y (or A) becomes infinitely large. 

For y4(l) = — 1, as ?/ becomes large / decreases and tends to a constant value, while A 
beginning at a negative value approaches zero eventually, tend to zero value as ?/ — oo. 
The Ricci scalar R for both initial conditions on A is non-singular at all finite y and decays 
with y monotonically. 

It is necessary to look into the behaviour of the quantity a^R as we change the initial 
values of /, A and r^. Firstly, let us note that changing the initial value of / does not 
produce any effect. However, interesting features do arise when we change A{\) or rc(l). 
FigS] shows this variation for different initial values of A. Starting from the top left, we 
change ^4(1) from a small value to successive larger values, while keeping rc(l) and /(I) 
fixed. We note that the location of the maximum shifts towards vertical axis and the value 
of a^R at a given y becomes larger with increasing initial A. Also, it is worth noting that 
the a^R is always positive and goes to zero asymptotically. 

In Figini we show how a variation of rc(l) (at fixed /(I) and ^4(1)) changes the behaviour 
of (j^R. At smaller values of rc(l) the flow continues indefinitely into the future. But for a 
critical value of rc(l), the flow develops a future singularity. Moreover, with changing 
we switch from positive to negative values of a^R. The flow cannot be extended beyond a 
certain value of y if the initial is above a certain critical value. 

To understand the appearance of genuine singularities, we investigate when such future 
singularities arise in the flow and their dependence on the initial conditions on r^ and A,. 
We search the range of 1 < y < 100 with different initial conditions on A and r^, and 
look for the appearance of singularities. Since, the flow does not depend on /, we can 
make a phase diagram in the space of (/1(1), rc(l)), in which we can represent singular and 
non-singular flows. Fig l6] shows such a plot. The white (blank) region represents conditions 
that do not lead to singularities, whereas the darker shades of gray denote flows with value 
oi y = ys at which the flow becomes singular. Thus, for future singularities, increasing 
rc(l) makes yg smaller, and hence the flow can only be continued to shorter times in the 
future. These trends can also be seen in FigsJH and [51 Note that for A{1) < the flow has 
no singularities at all and exists for all time. Also, from the phase diagram in Figl6]we see 
that, there are regions of initial conditions A{l),rc{l) corresponding to both types of flow 
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— non-singular and only future singular. 



Given the above results, let us now try to understand the consequences a little more. 
From the solutions for /, A, and the expression for a^R (all of which are functions of 
y) we note that for fixed y = yo, all these quantities have a fixed value, i.e. f{yo) = fo, 
A{yo) = Aq, r^iyo) = r^o and R{yo) = cr^Riyo) = Ro- However, we also have y^ = ^ 
or a = Hence, for fixed y, the flow along A is equivalent to moving along the extra 
dimension. Changing a (or A) keeps the obtained /, A, and a^R unchanged at fixed y. 
However, R will change with changing a and at different a the value of the Ricci scalar will 
obviously be different. Thus, the solutions obtained represents Ricci flow along the extra 
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(a)future singularities 

FIG. 6: phase diagram showing non-singular and singular flows in space of (^(1), ?'c(l)) 

dimensional coordinate (proportional to \/~X at fixed y), such that the value of /, A and r^, 
as we move from one location to another in the extra dimension remains unaltered, though 
the value of the five dimensional Ricci scalar changes. 
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IV. REMARKS AND CONCLUSIONS 



We conclude with a summary and some remarks. 

We have looked into Ricci flow of unwarped and warped product manifolds. In the 
unwarped case we have focused on x S™, H" x H™ and S" x H™. Distinguishing features of 
the Ricci flow of such manifolds have been highlighted with reference to flxed points/curves, 
singularity formation, intersections which imply isotropisation at some A and dependence 
on initial conditions. We have also discussed multiple products and provided illustrative 
examples for the case x S" x S™. 

Further, inspired by warped braneworlds, we have looked into Ricci flow of such manifolds 
with arbitrary functional forms of the warp factor and the scale of the extra dimensions. We 
discuss two separate cases-one where these functions are separable (this being analytically 
solvable) and the other, the non-separable case where we use numerical methods to flnd the 
evolution. In the former case, we are able to flnd the line element exactly and also compute 
its curvature. We note that the spacetime is necessarily AdS for all A and this A can, in 
some way, be associated with the k (AdS curvature scale) in the warped braneworld scenario. 
Thus Ricci flow seems to lead to a RS type warp factor and the curvature of the spacetime 
remains AdS during evolution. In the more general case, we have analysed numerically the 
particular case of scaling solutions. We flnd that, for certain initial conditions, the evolution 
is regular throughout. We have shown that Ricci flow in A turns out to be equivalent to a 
flow along cr (at flxed y) and the /, A,rc and a^R at that value of y do not change with 
varying A (or cr), though the value of R does evolve. For other initial conditions, we note that 
genuine singularities may appear at flnite y. We also observe a transition from non-singular 
to singular solutions as we tune the initial conditions. 

Finally, we mention that we hope to look into more complex and realistic scenarios in the 
warped context, in future. In the unwarped case, we would like to arrive at more general 
conclusions for arbitrary multiple products by a more detailed analysis of the associated 
dynamical system. 
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